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February 2, 2008 

Abstract. High-precision analyses of supersymmetry parameters aim at reconstructing the funda- 
mental supersymmetric theory and its breaking mechanism. A well defined theoretical framework 
is needed when higher-order corrections are included. We propose such a scheme, Supersymmetry 
Parameter Analysis SPA, based on a consistent set of conventions and input parameters. A repos- 
itory for computer programs is provided which connect parameters in different schemes and relate 
the Lagrangian parameters to physical observables at LHC and high energy e + e~ linear collider 
experiments, i.e., masses, mixings, decay widths and production cross sections for supersymmetric 
particles. In addition, programs for calculating high-precision low energy observables, the density 
of cold dark matter (CDM) in the universe as well as the cross sections for CDM search exper- 
iments are included. The SPA scheme still requires extended efforts on both the theoretical and 
experimental side before data can be evaluated in the future at the level of the desired precision. 
We take here an initial step of testing the SPA scheme by applying the techniques involved to a 
specific supersymmetry reference point. 
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1 INTRODUCTION 

At future colliders, experiments can be performed in 
the supersymmetric particle sector ^|z]|2G], if realized 
in Nature, with very high precision. While the Large 
Hadron Collider LHC can provide us with a set of well- 
determined observables (HE], in particular masses of 
colored particles and precise mass differences of var- 
ious particle combinations, experiments at the Inter- 



national e + e~ Linear Collider ILC jZflBK] offer high- 
precision determination of the non-colored supersym- 
metry sector. Combining the information from LHC on 
the generally heavy colored particles with the informa- 
tion from ILC on the generally lighter non-colored par- 
ticle sector (and later from the Compact Linear Col- 
lider CLIC on heavier states) will generate a com- 
prehensive high-precision picture of supersymmetry at 
the TeV scale JT] . Such an analysis can be performed 
independently of specific model assumptions and for 
any supersymmetric scenario that can be tested in lab- 
oratory experiments. It may subsequently serve as a 
solid base for the reconstruction of the fundamental su- 
persymmetric theory at a high scale, potentially close 
to the Planck scale, and for the analysis of the micro- 
scopic mechanism of supersymmetry breaking |12II13| . 

The analyses will be based on experimental accura- 
cies expected at the percent down to the per-mil level 
[5lll4| . These experimental accuracies must be matched 
on the theoretical side. This demands a well-defined 
framework for the calculational schemes in perturba- 
tion theory as well as for the input parameters. The 
proposed Supersymmetry Parameter Analysis Conven- 
tion (SPA) [Sect. 2] provides a clear base for calculating 
masses, mixings, decay widths and production cross 
sections. They will serve to extract the fundamental 
supersymmetric Lagrangian parameters and the super- 
symmetry-brcaking parameters from future data. In 
addition, the renormalization group techniques must 
be developed for all the scenarios to determine the 
high-scale parameters of the supersymmetric theory 
and its microscopic breaking mechanism. 

By constructing such a coherent and unified basis, 
the comparison between results from different calcula- 
tions can be streamlined, eliminating ambiguous pro- 
cedures and reducing confusion to a minimum when 
cross-checking results. 

A program repository [Sect. 3] has therefore been 
built in which a series of programs has been collected 
that will be expanded continuously in the future. The 
programs relate parameters defined in different schemes 
with each other, e.g. pole masses with DR masses, and 
they calculate decay widths and cross sections from the 
basic Lagrangian parameters. An additional set of pro- 
grams predicts the values of high-precision low-energy 
observables of Standard Model (SM) particles in su- 
persymmetric theories. The program repository also 
includes global fit programs by which the entire set 
of Lagrangian parameters, incorporating higher-order 
corrections, can be extracted from the experimental 



observables. In addition, the solutions of the renormal- 
ization group equations are included by which extrapo- 
lations from the laboratory energies to the Grand Uni- 
fication (GUT) and Planck scales can be performed 
and vice versa. Another category contains programs 
which relate the supersymmetry (SUSY) parameters 
with the predictions of cold dark matter in the uni- 
verse and the corresponding cross sections for search 
experiments of cold dark matter (CDM) particles. 

It is strongly recommended that the programs avail- 
able in the repository adopt the structure of Ref. \To\ 
for the Lagrangian, including flavor mixing and CP 
phases, and follow the generally accepted Supersym- 
metry Les Houches Accord, SLHA, for communication 
between different programs |16) . For definiteness, we 
reproduce from the supcrpotcntial (omitting R- 
parity violating terms), in terms of supcrficlds, 



W 



tab 



(Ys^H^Ej + {Y D ) ij H a d Q\D j 

h (Yu^HtQtUj - nHSHt] , (1) 



where the chiral superfields of the Minimal Supersym- 
metric Standard Model (MSSM) have the following 
SU{3) C <S> SU(2) L ® U(1) Y quantum numbers 



L: (1,2, -1), E: (1,1,1), Q: (3,2, ±), U : (3, 1,-f) 
D : (3, 1,|), H d : (l,2,-±), H u : (1,2,±). 

The indices of the SU(2)l fundamental representation 
are denoted by a, b = 1,2 and the generation indices by 
i,j = 1,2,3. Color indices are everywhere suppressed, 
since only trivial contractions are involved. e a b is the 
totally antisymmetric tensor, with 612 = e 12 = 1. 
The soft SUSY breaking part is written as 



Csoft = £ a b (T E )ijH%L iL e* jR + (T D ) i:j H d Q b lL d* R 



- (Tu ) ij Q\ h Uj R 



h.c. 



m 2 H H* da H a d 



^KaK - {mle ab H a d H b u + h.c.) 

+ u lR {ml) i3 u* R + d lR {m 2 i ) i] d* R + e 2R (ml) zj e* R 

+ i (Mxlb + M 2 w A w A + M 3 g x g x ^j + h.c. , (2) 



where the Hi are the scalar Higgs fields, the fields with 
a tilde are the scalar components of the supcrficld with 
the identical capital letter; the bino is denoted as b, 
the unbroken SU(2)l gauginos as w A=1 ' 2 ' 3 , and the 
gluinos as g x=1 " 8 , in 2-component notation. The T 
matrices will be decomposed as Ty = AijYij, where Y 
are the Yukawa matrices and A the soft supersymmetry 
breaking trilincar couplings. 

Much work on both the theoretical and the exper- 
imental side is still needed before data could be eval- 
uated in the future at the desired level of accuracy. 
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SPA CONVENTION 

— The masses of the SUSY particles and Higgs bosons are defined as pole masses. 

— All SUSY Lagrangian parameters, mass parameters and couplings, including tan j3, are given 
in the DR scheme and defined at the scale M = 1 TeV. 

— Gaugino/higgsino and scalar mass matrices, rotation matrices and the corresponding angles 
are defined in the DR scheme at M, except for the Higgs system in which the mixing matrix is 
defined in the on-shell scheme, the momentum scale chosen as the light Higgs mass. 

— The Standard Model input parameters of the gauge sector are chosen as Gf, a, Mz and 
a^ s {Mz)- All lepton masses are defined on-shell. The t quark mass is defined on-shell; the 
b, c quark masses are introduced in MS at the scale of the masses themselves while taken at a 
renormalization scale of 2 GeV for the light u, d, s quarks. 

— Decay widths/branching ratios and production cross sections are calculated for the set of pa- 
rameters specified above. 



Table 1. Definition of the supersymmetry parameter convention SPA 



These tasks of the SPA Project will be defined in de- 
tail in Sect. 4. 

In Sect. 5 we introduce the SUSY reference point 
SPSla' as a general setup for testing these tools in 
practice. This reference point is defined at a charac- 
teristic scale of 1 TeV in the Minimal Supersymmetric 
Standard Model with roots in minimal supergravity 
(mSUGRA). The point is a derivative of the Snow- 
mass point SPSla |17j : its parameters are identical 
except for a small shift of the scalar mass parameter 
and a change of the trilinear coupling to comply with 
the measured dark matter density ^H- Note, that the 
SPSla' parameters are compatible with all the avail- 
able high- and low-energy data. The parameters are 
close to point B' of Ref. \H)\. The masses are fairly 
light so that stringent tests of all aspects in the pro- 
gram can be performed for LHC and ILC experiments. 
The final target arc predictions on the accuracies of 
the fundamental supersymmetry parameters that can 
be expected from a common set of information when 
LHC and ILC experiments arc analyzed coherently. 

Additional benchmark points within and beyond 
mSUGRA, representing characteristics of different sce- 
narios, should complement the specific choice of SPSla'. 

2 SPA CONVENTION 

Extending the experience collected in analyzing Stan- 
dard Model parameters at the former e + e~ colliders 
LEP and SLC, we propose the set of conventions de- 
fined in Table ^ These conventions conform with the 
general SLHA scheme \Ijg\ but they are more specific in 
several points. 

Though largely accepted as standard, some of the 
definitions proposed in this SPA Convention should be 
explained in a few comments. 

For the SUSY Lagrangian parameters the DR sche- 
me |20II21) is most useful. It is based on regularization 
by dimensional reduction together with modified min- 
imal subtraction. This scheme is designed to preserve 



supersymmetry by maintaining the number of degrees 
of freedom of all fields in D dimensions, and it is tech- 
nically very convenient. The (3- functions for SUSY pa- 
rameters in this scheme are known up to 3-loop order 
|22| . It has recently been shown [23] that inconsisten- 
cies of the original scheme [21] can be overcome and 
that the DR scheme can be formulated in a mathe- 
matically consistent way. The ambiguities associated 
with the treatment of the Levi-Civita tensor can be 
parameterized as renormalization scheme dependence 
as was argued in |25| . Checks by explicit evaluation 
of the supersymmetric Slavnov- Taylor identities at the 
one-loop level have shown that the DR method gen- 
erates the correct counter terms [2H]. [We will use the 
version of the DR scheme as given in there re- 
ferred to as DR scheme.] To make use of the highly 
developed infrastructure for proton colliders, which is 
based on the MS factorization scheme [2Z|, a dictionary 
is given in Sect. 3. 2 for the translation between the DR 
and MS schemes, as well as the on-shell renormaliza- 
tion schemes. 

The SUSY scale is chosen M = 1 TeV to avoid 
large threshold corrections in running the mass pa- 
rameters by renormalization group techniques from the 
high scale down to the low scale. Fixing the scale M 
independent of parameters within the supersymmetry 
scenarios is preferable over choices relating to specific 
parameters, such as squark masses, that can be fixed 
only at the very end. By definition, this point can 
also be used to characterize uniquely multiple-scale ap- 
proaches. 

Mixing parameters, in particular tan/3, could have 
been introduced in different ways |29| : however, choos- 
ing the DR definitions proposed above has proven very 
convenient in practical calculations. 

The masses of Higgs bosons [3D] , in the MSSM of 
the charged H , of the neutral CP-odd A, and of the 
two CP-even h, H particles, are understood as pole 
masses, M H ± A Hh . For given Ma, the pole masses 
Mff,h of the CP-even Higgs bosons arc obtained as 
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poles q 2 = Mfj h of the dressed propagator matrix, 

a f 2^ = fq 2 -m 2 H + S HH (q 2 ) Z hH {q 2 ) \ 
Hh[q ' \ S hH (q 2 ) q 2 -m 2 +S hh (q 2 )) 

involving the tree-level masses mH,h and the diagonal 
and non-diagonal on-shell-renormalized self-energies E. 
In the on-shcll scheme, the input parameters are renor- 
malized on-shcll quantities, in particular the A-boson 
mass, with accordingly defined counter terms. 

Owing to the momentum dependence of the self- 
energies, there is no unique mixing angle (a) for the 
neutral CP-even Higgs system beyond the tree level, 
and the SPA choice can be understood as a convention 
for an "improved Born approximation" . A convenient 
choice for q 2 in the self-energies which minimizes the 
difference of such an approximation with respect to cal- 
culations involving the proper self-energies in physical 
matrix elements, is given by q 2 = M 2 . 

The physical on-shcll masses arc introduced in the 
decay widths and production cross sections such that 
the phase space is treated in the observables closest to 
experimental on-shell kinematics. This applies to the 
heavy particles while the masses of the light particles 
can generally be neglected in high energy processes. 

In the chargino/ncutralino sector the number of ob- 
servable masses exceeds the number of free parameters 
in the system, gaugino/higgsino mass parameters and 
tan (3. The most convenient set of input chargino/neu- 
tralino masses is dictated by experiment [the three low- 
est mass states in this sector, for example] while the 
additional masses are subsequently predicted uniquely. 
Similar procedures need to be followed in the sfermion 
sector. 

3 PROGRAM BASE 

3.1 PROGRAM CATEGORIES 

The computational tasks that are involved in the SPA 
Project can be broken down to several categories. Each 
of the codes that will be collected in the SPA program 
repository is included in one or more of these cate- 
gories. It is understood that in each case the theoretical 
state-of-the-art precision is implemented. For commu- 
nication between codes SLHA ^B] is strongly recom- 
mended, which is extended in a suitable way where 
appropriate. 

1) Scheme translation tools: 

The communication between codes that employ dif- 
ferent calculational schemes requires a set of trans- 
lation rules. In the SPA program repository we there- 
fore collect tools that implement, in particular, the 
definitions and relations between on-shell, DR and 
MS parameters in the Lagrangian as listed in Sect. 
3.2 below. 

2) Spectrum calculators: 

This category includes codes of the transition from 
the Lagrangian parameters to a basis of physical 



particle masses and the related mixing matrices. 
This task mainly consists of deriving the on-shell 
particle masses (including higher-order corrections) 
and of diagonalizing the mixing matrices in a con- 
sistent scheme, making use of the abovementioned 
tools as needed. 

3) Calculation of other observables: 
3A) Decay tables: 

compute the experimentally measurable widths 
and branching fractions. 

3B) Cross sections: 

calculate SUSY cross sections and distributions 
for LHC and ILC. 

3C) Low-energy observables: 

compute the values of those low-energy, high- 
precision observables [e.g., b — ► s'j, B s — > /i/i, 
gy, — 2] that are sensitive to SUSY effects. 

3D) Cosmological and astrophysical aspects: 

this category of programs covers the derivation 
of cold dark matter (CDM) relic density in the 
universe, cross sections for CDM particle search- 
es, astrophysical cross sections, etc. in the SUSY 
context. 

4) Event generators: 

Programs that generate event samples for SUSY 
and background processes in realistic collider envi- 
ronments. 

5) Analysis programs: 

These codes make use of some or all of the above 
to extract the Lagrangian parameters from experi- 
mental data by means of global analyses. 

6) RGE programs: 

By solving the rcnormalization-group equations, the 
programs connect the values of the parameters of 
the low-energy effective Lagrangian to those at the 
high-scale where the model is supposed to match to 
a more fundamental theory. High-scale constraints 
arc implemented on the basis of well-defined theo- 
retical assumptions: gauge coupling unification, 
mSUGRA, GMSB, AMSB scenarios, etc. 

7) Auxiliary programs and libraries: 

Structure functions, bcamstrahlung, numerical 
methods, SM backgrounds, etc. 

This is an open system and the responsibility for 
all these programs remains with the authors. SPA pro- 
vides the translation tables and the links to the com- 
puter codes on the web-page 

http : / / spa . desy . de/ spa/ 

Conveners responsible for specific tasks of the SPA 
Project will be listed on this web-page; the informa- 
tion will be routinely updated to reflect the momentary 
state of the project at any time. 

3.2 SCHEME TRANSLATION 

This subsection presents a few characteristic examples 
of relations between on-shell observables and DR, MS 
quantities at the electroweak scale Mz and the SUSY 
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scale M. For brevity, here only the approximate one- 
loop results are given PJ; it is understood that the 
codes in the program repository include the most up- 
to-date higher-loop results. 

(a) Couplings: 

• gauge couplings: 



ms _ dr , _ m_j_ r 
91 1 96^ ° ! 



(3) 



• Yukawa couplings between the gaugino A, ; the chi- 
ral fermion ipk and the scalar </>&: 



-MS _ 

9ik — 9i 



DR I , . (ffP R ) n \ "* \3l ) n i 



(g? B f 



(4) 



• Yukawa couplings between the scalar 4>i and the 
two chiral fermions ipj and ipk'- 

Y™ = Y™U + ^J^l [Cf'-2CF+Cp]) (5) 
^ i=i ' 

• trilinear scalar couplings: 

These couplings do not differ in the two schemes. 

Ci and CI are the quadratic Casimir invariants of 
the adjoint representation and the matter represen- 
tation r of the gauge group Gi, respectively. They 
are given by C t = [3,2,0] for [St/(3), SU(2), U(l)} 
and CT = [4/3,3/4,3/5 x Y r 2 } for the fundamental 
representations of SU(3), SU(2), and the U(l) hy- 
percharge Y r . 

(b) SUSY PR, MS and pole masses: 

• gaugino mass parameters 



M 



MS 



At 



16tt 2 



•a 



(6) 



• higgsino mass parameter. 

' 2 




(7) 



C, ff denoting the SU(2) and f7(l) Casimir invari- 
ants of the Higgs fields. 



where S denotes the fermion self-energy renormal- 
ized according to the DR-scheme at the scale M. 
As an explicit example we note the one-loop re- 
lation between the SU(3) gaugino mass parameter 

M3(M) DR and the gluino pole mass m§ [without 
sfcrmion mixing] at the one-loop order: 



irig = M 3 



DR (M) 

47T 



(9) 



<K 



m g [ 15 + 9 In 



M2 
2 



q 2—1 



m g B 1 (m 2 g ,mg, 



where Bi is the finite part of one of the one-loop 
two-point functions at the scale in the DR scheme 
M (and analogously A ,B to be used later), cf. 

Ref. Eg. 

• scalar pole masses: 

A similar relation holds for the squared scalar masses 



= M 2 ' BR 



m i, pole 



S(q 2 = ml polc ) (10) 



The one-loop QCD corrections for the left squarks 
of the first two generations in the limit of vanishing 
quark masses may serve as a simple example: 



m§ = A£ DR (Af) 



(11) 



2a" H (M) 
3tt 



(m? - rn~)B (rr4,rn~,0) 



2m?£ (m?, rof, 0) + A (?n?) - A (m?) 



(c) SM parameters: 

The following paragraphs summarize the SM input 
values for the analysis. Only approximate formulae 
are presented for brevity, while the complete set of 
relations is available on the program repository. 
In a few cases the evolution from the scale Mz 
to M is carried out by means of RGEs instead of 
fixed-order perturbation theory because they have 
proven, presently, more accurate; this may change 
once the necessary multi-loop calculations will be 
completed. 



• sfermion mass parameters: 

These parameters do not differ in the DR and MS 
schemes. 



• fermion pole masses: 

The pole masses can be written schematically as 
mj, pole = M® R - Re E = m h polo ) (8) 



a, a (Mz), <#${M) 
a m (M z ) 



1 - Aa SM - A^susy 



(12) 



In 



4£ln 



6tt *~ M z ' * M z 

2—1 

2 

N,.Ol In 

M z 



/ »=i 
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Aa$M summarizes the SM contributions from the 
lcptons, quarks and the VU-boson. In the SUSY 
contributions, Z\asuSY, / sums over all charged 
sfermions, N c is the color factor and Q / the (s)ferm- 
ion charge. 



„DR 



(M z ) 



, DR 



(M) 



! 9 DR {M Z 



1 



„DR 



(M) 



^ R {M Z ) ^Mf 
47rcos 2 6» DR (Afz) M 2 

a m jAIz) 
sin 2 e m {M z ) 



(13) 



1 



„DR 



(M z 



_ ln M| 
4?r S in 2 9 m {M z ) n M 2 



(14) 



• sin 2 8 UR at M z and at M: 



The electroweak mixing parameter sin 6 R (Mz) is 



given by 

sin 2 OM z ) 



1 



sin a {AJ-z, 

vDR, 



ira 



l (Mz) 



y/2M 2 G F {\ - Af) 



(15) 



where the contributions from loops of SM and SUSY 
particles are denoted by Af |33U34| . At the scale M 
the electroweak mixing parameter can be calculated 
subsequently from 

tan 2 e m {M) = a™ {M) /af*(M) (16) 

by making use of the couplings af R (M) given in 
the preceeding paragraph. 



,2 /)DR 



and sin 8 c fi at M 



The electroweak mixing angle in the effective lep- 
tonic (electronic) vertex of the Z boson is defined 
as 



SB in a *W(Jlfc) = i(l-Reg] 



in terms of the effective vector and axial vector cou- 
plings g v A of the Z to electrons. The relation to 

sin 2 9 DR (Mz) is given by (at one-loop order) 



sin 2 9 DR {M Z ) = sin 2 6» e ff 

n lZ {M 2 ) 



sin 2 6. 



eg 



il 7 z(0) 



2 Ml 



(18) 



involving the photon- Z* non-diagonal self-energy 
n^zi.0 2 ) and the non- universal electron- Z vertex 
correction form factors f v A (q 2 ), 

/ c = \ fv(M z ) - a - 2 sin 2 9 cS ) f e A (M 2 z ), (19) 



with all the loop quantities renormalized in the DR 
scheme at the scale Mz- For explicit expressions 
see |53I51] . 



• a 



DR 



ai M z and M, related to af h (M z ) 



, DR 



„MS 



(M z ) 



Aa, = 



(M z 



l-Aa s 
a s (M z ) 



(20) 



2tt 



1 _ 2 ln TO * 

2 3 n M z 



2m ^ _ iVVln^i 
A / ^ Z^ Z^ 



DR 



(M) 



M z 
vF(Mz) 



£a?*(M z ) In §f 



(21) 



• W, Z bosons, pole and DR masses: 



The pole masses M v (V = W 7 Z) and the DR 
masses at Mz are related by 



M, 



M 2 ' DR (M Z ) 



Ren^ v (p 2 = M 2 ) (22) 



involving the renormalized transverse vector-boson 
self-energies in the DR scheme at the scale Mz- The 
Z pole mass is a direct input parameter, whereas 
the W pole mass is derived from the relation to the 
low-energy parameters a and Fermi constant Gf 
according to the SPA Convention: 



Mw 



1 - 



M 



M 2 



V2G F (1 - At 



(23) 



Ar summarizes the loop contributions from the SM 
and SUSY particles as given explicitly in 



The self-energies at the scale M can be written 
symbolically as 



16tt 2 7J| z = 16tt 2 7T| Zj SM +Hi g 



(24) 



E 



mfv 2 z, l] B 22 (M 2 ,m 2 h ,m)) 



[fijzH(M%, , ) 



+ 2g ijZ B (M%, 77%, m^)] 



i6ir 2 nT 



WW 



\Q-K 2 II WW ^ SM+Higgs 



(25) 



- Y, - V . f '1u,, B 22 (M w , mj , mj, ) 



+ 2g i jwBo(Mw,m jt o,m~+)] 

where Vfvij are the couplings of the gauge boson 
to sfermions and fyv and gijv are combinations of 
left- and right-couplings to charginos and neutrali- 
nos; B22 and H are combinations of the Bi and A4 
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loop functions. Detailed formulae are given in |36| . 

• charm and bottom running MS mass at m cb and 
DR mass at M z , cf. [37H3ffl 



,MS 



(m b ) 



a 



MS 



(M 2 



a 

DR 



37T 



* s (ro 6 ) 

23a 2 - DR 
72tt 



mf K (M z ) 



1 - Am b (M z ) 



(26) 



(27) 



2a, 



Am h (M z ) = ^rrigfj. tan /3 /(m^ , m~ a , to?) 



If 



/(« 2 ,& 2 ,c 2 ) 



+ ^-2^tMtan/3/(m fi ,m t - 2 ,/i ) 

2 

- Y^^2Mtan/3 

x [ cos 2 6~ t I (ml , M| , M 2 ) + I {f^fo} 
+ {cos — > sin; Qi — ► Q2}] 

a 2 6 2 loga 2 /& 2 + cyclic 
(a 2 - 6 2 )(6 2 - c 2 )(a 2 - c 2 ) 



with X£(M Z ) = 27 6 (Afz) - m^ R (M z )Am b (M z ) 
and E b (M z ) being the self-energy of the bottom 
quark due to supcrsymmetric particles and heavy 
SM particles and Am b (M z ) including the large fi- 
nite terms proportional to tan (3 which have been 
resummed |38| . In the case of the charm quark the 
additional running between m c and m b has to be 
included. The SUSY contributions are in general 
small and no resummation is necessary. The masses 
are evolved from the scale M z to M by means of the 
RGEs for the Yukawa couplings as described below. 

• top quark pole mass and DR mass at M z : 



DR. 



{M z ) = ro t 



5a 



DR 



„DR. 



-Ci 



3tt 



log 



- E 



(28) 



where ct(M\ / 'to 2 ) is the gluonic two-loop contribu- 
tion and E accounts for the electroweak as well as 
the SUSY contributions. The mass is evolved to the 
scale M by means of the Yukawa RGEs; see next. 

• Yukawa couplings and running masses of SM par- 
ticles at M : 

The vacuum expectation values v® R and are 
initially given by: 



M^(M Z ) 



- t 9 



2 : DR 



[M Z 



(29) 



2.DR 



(M Z ) 



2, DR. 



(M z ) 



v» R (M z )/v% R (M z ) = tan /? DR (Mz) (30) 

tan/3 DR (Mz) must be evolved down from the con- 
ventional parameter tan /? DR (M) by means of RGE. 
From the DR masses at M z the Yukawa couplings 
are calculated: 

Y™(M Z ) = V2m™(M z )/v™(M z ) (31) 



Y b U T R (M z ) 



V2m? R (M z )/v™(M z ) (32) 



In a second step, they are evolved together with 
the gauge couplings and the vacuum expectation 
values to M via RGEs. At this scale the running 
SM fermion masses and gauge boson masses are re- 
lated to the Lagrangian parameters by the usual 
tree-level relations. This is, presently, a better ap- 
proach for the evolution of the Yukawa couplings 
than fixed-order perturbation theory. 

3.3 WIDTHS AND CROSS SECTIONS 

(a) Decay widths: 

The decay widths are defined as inclusive quanti- 
ties including all radiative corrections; the masses 
of the heavy particles are taken on-shcll, light par- 
ticle masses are set zero. 

(b) Cross sections for e + e~ collisions: 

Cross sections, cr(e + e~ — > {F}), for the production 
of a set of supersymmetric particles/Higgs bosons 
{F} are defined at the experimental level in e + e~ 
collisions including up-to-date radiative corrections 
except hard 7 bremsstrahlung to exclude large con- 
tributions from radiative return. 
In general, large QED-typc photonic corrections 
cannot be disentangled from genuine SUSY-specific 
parts, and in the comparison of theoretical predic- 
tions with experimental data all higher-order terms 
have to be included. To elucidate the role of the spe- 
cific supcrsymmetric loop corrections, a reasonable 
and consistent prescription for cut-independent re- 
duced cross sections shall therefore be defined. Since 
the leading QED terms arising from virtual and real 
photon contributions that contain large logarithms 
can be identified and isolated, the "reduced" gen- 
uine SUSY cross sections are defined, at the theo- 
retical level, by subtracting the logarithmic terms 
log4/ii5 2 /s in the soft-photon energy cut-off AE 
and in log s / m 2 from non-collinear and collinear 
soft 7 radiation off light fcrmions / = e, /i, . . . and 
virtual QED corrections. In this definition of re- 
duced cross sections [see also the logarithmi- 
cally large QED radiative corrections are consis- 
tently eliminated in a gauge-invariant way. By the 
same token, the reduced cross sections are defined 
without taking into account bcamstrahlung. 
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(c) Cross sections for hadron collisions: finite shift from the commonly used MS density func- 

Cross sections for proton collisions at Tevatron and tions to thc DR density functions has to be carried 

LHC, cr(pp {F}), include all QCD and other out 02]. Moreover, for massive final state particles spu- 

available corrections, with infrared and collinear rious density functions for the (4 - D) gluon compo- 

singularities tamed by defining inclusive observables, nents have to be introduced to comply with the factor- 

or properly defined jet characteristics, and intro- ization theorem, see 03E2] for details. Formulating an 

ducing the renormalized parton densities, provided efficient combination of the most attractive elements of 

parametrically by the PDF collaborations gOETj. both schemes in describing hadronic processes is there- 
fore an important task of the project. 



4 TASKS OF THE SPA PROJECT 

A successful reconstruction of the fundamental struc- 
ture of the supersymmetric theory at the high scale 
and the proper understanding of the nature of cold 
dark matter from experimental data require thc precise 
analysis of all information that will become available 
from collider experiments, low-energy experiments, as- 
trophysical and cosmological observations. Preliminary 
studies [see Sect. 5], initiating this SPA Project, have 
shown that while this aim can in principle be achieved, 
it still needs much additional work both on the theoret- 
ical as well as on the experimental side. In particular, 
we identify the following areas of research as central 
tasks of the SPA Project: 

Higher-order calculations 

While the precision of SUSY calculations has gradually 
shifted from leading-order (LO) to next-to-leading or- 
der (NLO) accuracy [and, in some areas, beyond], the 
present level still does not match the expected exper- 
imental precision, particularly in coherent LHC+ILC 
analyses. Thc experimental precision, however, has to 
be fully exploited in order to draw firm conclusions on 
the fundamental theory. To close this gap, the SPA 
Project foresees new efforts to push the frontier in 
higher-order SUSY calculations to the line necessary 
for the proper interpretation of experimental analyses. 

Improving the understanding of thc DR scheme 

The DR scheme recommended for higher-order calcula- 
tions can be formulated in a mathematically consistent 
way \2'6\ and is technically most convenient. Many ex- 
plicit checks at the one-loop level have shown that the 
DR method generates the correct counter terms. How- 
ever, there is no complete proof yet that it preserves su- 
persymmetry and gauge invariance in all cases. There- 
fore, as the precision of SUSY calculations is pushed 
to higher orders, the SPA Project also requires fur- 
ther investigation of the symmetry identities in thc DR 
scheme. 

Moreover, there is an obvious dichotomy between 
the DR scheme, which is convenient for the definition of 
SUSY parameters and their rcnormalization group evo- 
lution, and the MS scheme, which is generally adopted 
for thc calculation of hadronic processes [27j. While, as 
argued before, the MS scheme requires ad-hoc counter 
terms to restore supcrsymmctry, in the DR scheme a 



Improving experimental and theoretical precision 

The set of observables that has been included so far 
in experimental analyses, by no means exhausts the 
opportunities which data at LHC and at ILC are ex- 
pected to provide in the future. SPA Project studies 
will aim to identify any new channels that can give ad- 
ditional information, cither independent or redundant 
[improving fit results], and they will include them in 
a unified framework. In connection with realistic es- 
timates of theoretical uncertainties, a solid account 
of error sources and correlations has to be achieved. 
Furthermore, the sophistication of thc experimental 
results will be refined by including more precise sig- 
nal and background calculations, and improved simu- 
lations as mandatory for the analysis of real data. 

Coherent LHC + ILC analyses 

We put particular emphasis on the coherent combina- 
tion of future data obtained at LHC and ILC. While 
the LHC will most likely discover SUSY particles, if 
they exist, and will allow for the first tests of the SUSY 
paradigm, e + e~ data make possible high-precision in- 
vestigations of the weakly-interacting sector. Feedback 
and coherently combined analyses, which will greatly 
benefit from a concurrent running of both colliders, are 
indispensable for a meaningful answer to thc questions 
raised in thc present context. Studies as initiated by 
the LHC/LC Study Group 05] are a vital part of the 
SPA Project. 

Determining SUSY Lagrangian parameters 

While at leading order the Lagrangian parameters con- 
nected with different supersymmetric particle sectors 
can in general be isolated and extracted analytically 
from closely associated observables, the analysis is much 
more complex at higher orders. Higher orders intro- 
duce the interdependence of all sectors in the observ- 
ables. The development of consistent analyses for the 
global determination of the Lagrangian parameters in 
this complex situation has started and, conform with 
general expectations for iterative steps in perturbative 
expansions, they can be carried out consistently with 
as few assumptions as possible. The set of Lagrangian 
parameters and their experimental error matrix can be 
determined, including higher-order corrections. How- 
ever, the experimental procedure must still be supple- 
mented by corresponding theoretical errors and their 
correlations. 



8 



Supersymmetry Parameter Analysis: SPA Convention and Project 



Cold dark matter 

As the precision is refined, astrophysical data play an 
increasingly important role in confronting supersym- 
metry with experiments. The class of models conserv- 
ing R-parity predict a weakly interacting, massive, sta- 
ble particle. The relic abundance of this particle im- 
poses crucial limits on supersymmetric scenarios 
While among the supersymmetry breaking models ver- 
sions of mSUGRA and of gaugino mediation 07J have 
been analyzed in detail, the analyses have to be ex- 
tended systematically to other scenarios. In models 
that account for the relic density, specific requirements 
on the accuracies must be achieved when the CDM 
particle is studied in high-energy physics laboratory 
experiments In turn, predictions based on the 

comprehensive parameter analysis of high-energy ex- 
periments determine the cross sections for astrophysi- 
cal scattering experiments by which the nature of the 
cold dark matter particles can be established. The SPA 
Project provides a platform for a systematic and con- 
tinuous interplay between the astrophysics and high- 
energy physics disciplines and the mutual refinement 
of their programs in the future. 

Extended SUSY scenarios 

The MSSM, in particular the parameter set SPSla' 
that we suggest for a first study, provides a benchmark 
scenario for developing and testing the tools needed for 
a successful analysis of future SUSY data. However, 
neither this specific point nor the MSSM itself may be 
the correct model for low-scale SUSY. Various param- 
eter sets [for instance other representative mSUGRA 
points as well as non-universal SUGRA, GMSB, AMSB, 
and other scenarios, see Ref. |49) for a brief summary] 
and extended models have therefore to be investigated 
within the SPA Project. In particular, models which 
incorporate the right-handed neutrino sector, must be 
analyzed extensively |50| . Furthermore, CP violation, 
i?-parity violation, flavor violation, NMSSM and ex- 
tended gauge groups are among the roads that nature 
may have taken in the SUSY sector. The SPA con- 
ventions are formulated so generally that they can be 
applied to all these scenarios. The goal of deriving the 
fundamental structure from data will also to be pur- 
sued for many facets in this more general context. 

5 EXAMPLE: REF POINT SPSla' 

To test the internal consistency of the SPA scheme and 
to explore the potential of such extended experimental 
and theoretical analyses we have defined, as an exam- 
ple, the CP and R-parity invariant MSSM reference 
point SPSla'. Of course, the SPA Convention is set up 
to cover also more general scenarios. 

The results for SPSla' presented below are based on 
preliminary experimental simulations. In some cases, 
however, extrapolations from earlier analyses for SPSla 
and other reference points have been used in order 
to substitute missing information necessary for a first 



Parameter 


SM input 


Parameter 


SM input 


m e 


5.110 • itr 4 


m^ oie 


172.7 


"V 


0.1057 


m b (m b ) 


4.2 


m T 


1.777 


m z 


91.1876 


m u (Q) 


3- 10~ 3 


G F 


1.1664 ■ 10~ 5 


m d (Q) 


7- 10 -3 


1/ct 


137.036 


m s (Q) 


0.12 




0.02769 


m c (m c ) 


1.2 




0.119 



Table 2. Numerical values of the SM input to SPSla' . 
Masses are given in GeV, for the leptons and the t quark 
the pole masses, for the lighter quarks the MS masses either 
at the mass scale itself, for c, b, or, for u, d, s, at the scale 
Q = 2 GeV. 



comprehensive test of all aspects of the SPA Project. It 
is obvious that many detailed simulations are needed 
to demonstrate the full power of predicting the funda- 
mental supersymmetric parameters from future sets of 
LHC and ILC data. 

In e + e~ annihilation experimental progress is ex- 
pected for the heavy chargino and neutralinos. Com- 
bining the results of such studies with LHC data ap- 
pear very promising and lead to improved mass deter- 
minations |51| . New techniques to determine slepton 
masses from cascade decays as very narrow resonances 
[521153] should be applied. For cross section measure- 
ments and other sparticle properties methods to deter- 
mine the decay branching ratios should be developed. 
At the LHC a recently proposed mass relation method 
offers substantial improvements in the reconstruction 
of squark and gluino masses |54j . 

Analysis of SUSY Lagrangian parameters 

The roots defining the Reference Point SPSla' are the 
mSUGRA parameters [in the conventional notation for 
CMSSM - see jSS] for the tighter original definition] in 
the set 



M 1/2 = 


250 GeV 


sign(^) = +1 


M = 


70 GeV 


tan/3(M) = 10 


= 


-300 GeV 





The left column, listing the universal gaugino mass 
M 1 / 2 , the scalar mass M and the trilinear coupling 
Aq [Yukawa couplings factored out], is defined at the 
GUT scale Mqut- The point is close to the original 
Snowmass point SPSla |17| : the scalar mass param- 
eter Mo is lowered slightly at the GUT scale from 
100 GeV to 70 GeV and A is changed from -100 GeV 
to —300 GeV. The values of the SM input parameters 
are collected in Table [5] Extrapolation of the above 
mSUGRA parameters down to the M = 1 TeV scale 
generates the MSSM Lagrangian parameters. Table 
displays the couplings and mass parameters after being 
evolved from Mqtjt to M using the RGE part of the 
program SPheno |5fc>j which is based on two-loop anal- 
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Parameter 


b-Pbla value 


Parameter 


ara la value 


9 


0.3636 


Mi 


103.3 


9 


0.6479 


M 2 


193.2 


9s 


1.0844 


M 3 


571.7 


Yr 


0.1034 


A T 


-445.2 


Y t 


0.8678 


A t 


-565.1 


Y b 


0.1354 


A b 


-943.4 




396.0 


tan f3 


10.0 


Mn d 


159.8 


\ M H U \ 


378.3 


M Ll 


181.0 


Ml 3 


179.3 


M El 


115.7 


M Ea 


110.0 


M Ql 


525.8 


Mq 3 


471.4 


M Vl 


507.2 


M U3 


387.5 


M Dl 


505.0 


M D3 


500.9 



Table 3. The DR SUSY Lagrangian parameters at the scale 
M = 1 TeV in SPSld from [mass unit in GeV; M% u 
negative]. In addition, gauge and Yukawa couplings at this 
scale are given in the DR scheme. 



Particle 


Mass [GeV] 


locale [GeV] 


h° 
H° 


116.0 
425.0 


1.3 
0.7 


-0 

Xi 

xl 

xi 
xf 


97.7 
183.9 
413.9 
183.7 


0.4 
1.2 
1.2 
1.3 


Sr 
eh 
fi 


125.3 
189.9 
107.9 


1.2 
0.4 
0.5 


h 
bi 


547.2 
564.7 
366.5 
506.3 


9.4 
10.2 
5.4 
8.0 


9 


607.1 


1.4 



Table 4. Supersymmetric masses for the SUSY scale M = 
1 TeV, and their variation if M is shifted to 0.1 TeV. 



yses of the /3-functions as well as the other evolution 
coefficients (other codes can be used equally well). 

This SPSla' set is compatible with all high-energy 
mass bounds and with the low-energy precision data, 
as well as with the observed CDM data, calculated as 
B(b -► s<y) = 3.0-10" 4 ^[ff-2] M /2 = 34-10" 10 EH, 
Aosusy = 2.1 ■ lCT 4 |58, and J?cdm/* 2 = 0.10 |SZj- 

The physical [pole] masses of the supersymmetric 
particles are presented in Table [S] The connection be- 
tween the Lagrangian parameters and the physical pole 
masses is presently encoded at the one-loop level for 
the masses of the SUSY particles, and at the two-loop 
level for the Higgs masses. QCD effects on the heavy 
quark masses are accounted for to two-loop accuracy. 



A systematic comparison with the other public pro- 
grams ISAJET jSnj, S0FTSUSY [HO] and SuSpect [HU has 

been performed in to estimate the technical accu- 
racy that can presently be reached in the evolution. 
The codes include full two-loop RGEs for all parame- 
ters as well as one-loop formulas for threshold correc- 
tions. The agreement between the actual versions of 
these calculations is in general within one percent. A 
special case are the on-shcll masses of the Higgs bosons 
which have been calculated by FeynHiggs |58| start- 
ing from the SPheno Lagrangian parameters as input. 
Here, discrepancies for the mass of the lightest Higgs 
boson amount to 2% or more which can be attributed 
to different renormalization schemes (see also jBH] for 
detailed discussions). 

Besides the comparison between different codes for 
spectrum calculations, a crude internal estimate of the 
theoretical errors at the present level of the loop calcu- 
lations may be obtained by shifting the matching point 
M from 1 TeV down to 0.1 TeV. A sample of parti- 
cle mass shifts associated with such a variation of the 
SUSY scale parameter is displayed in Tabled With er- 
rors at the percent level, the experimental precision at 
LHC can be matched in general. However, it is obvious 
that another order of magnitude, the per-mil level, is 
required in the theoretical precision to match the ex- 
pected experimental precision at ILC and in coherent 
LHC/ILC analyses - i.e., calculations of the next loop 
are called for 1 . 

To perform experimental simulations, the branch- 
ing ratios of the decay modes are crucial: these have 
been calculated using FeynHiggs [5H| and SDECAY j^Sl ; 
similar results may be obtained using CPSuperH |66| . 
The most important decay channels of the supersym- 
metric particles and Higgs bosons in SPSla' are col- 
lected in the Appendix, while the complete set is avail- 
able from the SPA web-site. Cross sections for the pro- 
duction of squarks, gluinos, gauginos and sleptons at 
the LHC are presented as a function of mass including 
the point SPSla'. Typical cross sections for pair pro- 
duction of charginos, neutralinos and sleptons at the 
ILC are presented for the point SPSla' function 
of the collider energy 

If SPSla', or a SUSY parameter set in the range of 
similar mass scales, is realized in nature, a plethora of 
interesting channels can be exploited to extract the ba- 
sic supersymmetry parameters when combining exper- 
imental information from sharp edges in mass distribu- 
tions at LHC with measurements of decay spectra and 
threshold excitation curves at an e + e~ collider with en- 
ergy up to 1 TeV |11| . From the simulated experimental 
errors the data analysis performed coherently for the 
two machines gives rise to a very precise picture of the 
supersymmetric particle spectrum as demonstrated in 
Table E] 

1 With p functions and evolution coefficients in the RGEs 
already available to third order |22l . the calculation of the 
two-loop order for the relation between the Lagrangian pa- 
rameters and the physical pole masses have been carried 
out in the approximation of massless vector bosons |64| 
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Particle 


Mass [GeV] 


Particle 


Mass [GeV] 


h° 
H° 
A 


116.0 
425.0 
424.9 


T\ 
T2 
V T 


107.9 
194.9 
170.5 


H+ 


432.7 


Ur 


547.2 


X? 

X.I 
X3 


97.7 
183.9 
400.5 


UL 

d R 
d L 


564.7 
546.9 
570.1 


xi 
xt 
xt 


413.9 
183.7 
415.4 


tl 

u 

bi 


366.5 
585.5 
506.3 


e~R 


125.3 


b 2 


545.7 


Sl 
i>e 


189.9 
172.5 


9 


607.1 




Table 5. Mass spectrum of supersymmetric particles 156V and Higgs bosons \58( in the reference point SPSla' . The 
masses in the second generation coincide with the first generation. 



Particle 


Mass 


"LHC" 


"ILC" 


"LHC+ILC" 


h° 


116.0 


0.25 


0.05 


0.05 


H° 


425.0 




1.5 


1.5 


X? 


97.7 


4.8 


0.05 


0.05 


xl 


183.9 


4.7 


1.2 


0.08 


xi 


413.9 


5.1 


3-5 


2.5 


xt 


183.7 




0.55 


0.55 


&R 


125.3 


4.8 


0.05 


0.05 


e L 


189.9 


5.0 


0.18 


0.18 




107.9 


5-8 


0.24 


0.24 




547.2 


7-12 




5-11 




564.7 


8.7 




4.9 


h 


366.5 




1.9 


1.9 


bi 


506.3 


7.5 




5.7 


9 


607.1 


8.0 




6.5 



Table 6. Accuracies for representative mass measurements 
of SUSY particles in individual LHC, ILC and coherent 
"LHC+ILC" analyses for the reference point SPSla' [mass 
units in GeV]. qR and <jl represent the flavors q = u, d, c, s. 
[Errors presently extrapolated from SPSla simulations.] 



While the picture so far had been based on evaluat- 
ing the experimental observables channel by channel, 
global analysis programs have become available 07 
EH] in which the whole set of data, masses, cross sec- 
tions, branching ratios, etc. is exploited coherently to 
extract the Lagrangian parameters in the optimal way 
after including the available radiative corrections for 
masses and cross sections. With increasing numbers of 
observables the analyses can be expanded and refined 
in a systematic way. The present quality of such an 



analysis |68j can be judged from the results shown in 
Tabled These errors are purely experimental and do 
not include the theoretical counterpart which must be 
improved considerably before matching the experimen- 
tal standards. 

Extrapolation to the GUT scale 

Based on the parameters extracted at the scale M, we 
can approach the reconstruction of the fundamental su- 
persymmetric theory and the related microscopic pic- 
ture of the mechanism breaking supersymmetry. The 
experimental information is exploited to the maximum 
extent possible in the bottom-up approach |12j in which 
the extrapolation from M to the GUT/Planck scale 
is performed by the renormalization group evolution 
for all parameters, with the GUT scale defined by the 
unification point of the two electroweak couplings. In 
this approach the calculation of loops and [3 functions 
governing the extrapolation to the high scale is based 
on nothing but experimentally measured parameters. 
Typical examples for the evolution of the gaugino and 
scalar mass parameters are presented in Fig. d While 
the determination of the high-scale parameters in the 
gaugino/higgsino sector, as well as in the non-colored 
slcpton sector, is very precise, the picture of the col- 
ored scalar and Higgs sectors is still coarse, and strong 
efforts should be made to refine it considerably. 

On the other hand, if the structure of the theory at 
the high scale was known a priori and merely the ex- 
perimental determination of the high-scale parameters 
were lacking, then the top-down approach would lead 
to a very precise parametric picture at the high scale. 
This is apparent from the fit of the mSUGRA parame- 
ters in SPSla' displayed in Table GO EH- A high-quality 
fit of the parameters is a necessary condition, of course, 
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Fig. 1. Running of the gaugino and scalar mass parameters as a function of the scale Q in SPSla' \5(tf . Only experimental 
errors are taken into account; theoretical errors are assumed to be reduced to the same size in the future. 



Parameter 


SPSla' value 


Fit error [exp] 


Mi 


103.3 


0.1 


M 2 


193.2 


0.1 


M 3 


571.7 


7.8 




396.0 


1.1 


M Ll 


181.0 


0.2 


Me 1 


115.7 


0.4 


Ml 3 


179.3 


1.2 


M Ql 


525.8 


5.2 


M Dl 


505.0 


17.3 


M Q3 


471.4 


4.9 


771A 


372.0 


0.8 


At 


-565.1 


24.6 


tan f3 


10.0 


0.3 



Table 7. Excerpt of extracted SUSY Lagrangian mass and 
Higgs parameters at the supersymmetry scale M = 1 TeV 
in the reference point SPSla! [mass units in GeV J. 



for the theory to be correct - however it is not a suffi- 
cient condition; deviations from the theory may hide in 
large errors of some observables which do not spoil the 
quality of the fit in the top-down approach but which 
are manifest in the bottom-up approach. 

Cold dark matter 

Constraints on SUSY cold dark matter can be obtained 
at LHC by specifying the underlying scenario and ana- 
lyzing all data simultaneously within the given bench- 
mark model. From a study of the SPSla point, based 
on very large statistics one may expect that the 
relic density can be determined to ~ 6% for the SPSla' 
scenario. For SPSla', the relic density depends on the 



Parameter 


SPSla' value 


Experimental error 


Mgut 


2.47 • 10 16 GeV 


0.02 ■ 10 ie GeV 


a GCT 


24.17 


0.06 


Mi 


250 GeV 


0.2 GeV 


Afo 


70 GeV 


0.2 GeV 


A 


-300 GeV 


13.0 GeV 


I 1 


396.0 GeV 


0.3 GeV 


tan/3 


10 


0.3 



Table 8. Comparison of the ideal parameters with the ex- 
perimental expectations in the top-down approach \68f . 



parameters of the neutralino and sfermion sector as the 
dominant channels are annihilation of neutralinos into 
fcrmion pairs and coannihilation with staus. In partic- 
ular, for the most sensitive component, coannihilation 
processes, the relic density is essentially given by the 
mass difference between the lightest slepton fi and the 
LSP Xii which can be directly measured at the ILC. 
Studies of fi production at threshold [23] and decay 
spectra to Xi m the continuum [7J suggest that for 
SPSla', even with moderate luminosity, a precision of 
~ 2% on the cold dark matter abundance is achievable. 
A systematic analysis of various scenarios is being car- 
ried out in the LCC project as well as by other 
groups. 



6 SUMMARY AND OUTLOOK 

If supersymmetry is realized in Nature, future experi- 
ments at the LHC and the ILC will provide very precise 
measurements of supcrsymmetric particle spectra and 
couplings. On the theoretical side these measurements 
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must be matched by equally precise theoretical calcu- 
lations and numerical analysis tools. The SPA Project, 
a joint theoretical and experimental effort, aims at pro- 
viding 



— a well-defined framework for SUSY calculations and 
data analyses, 

— all necessary theoretical and computational tools, 

— a testground scenario SPSla', 

— a platform for future extensions and developments. 



On this basis coherent analyses of experimental data 
can be performed and the fundamental supcrsymmct- 
ric Lagrangian parameters can be extracted. They can 
serve as a firm base for extrapolations to high scales 
so that the ultimate supersymmetric theory and the su- 
persymmetry breaking mechanism can be reconstructed 
from future data. 

Much work is still needed on the experimental and 
theoretical side to achieve these goals at the desired 
level of accuracy. Some of the short- and long-term 
subprojects have been identified and should be pursued 
in the near future. 

The SPA Project is a dynamical system expected 
to evolve continuously. The current status of the SPA 
Project, names of the conveners responsible for spe- 
cific tasks as well as links to the available calculational 
tools, can be found at the SPA home page 
|http : / / spa . desy . de/ spa/| 



APPENDIX 

(a) Decays of Higgs and SUSY particles in SPSla' 

The branching ratios of Higgs bosons and SUSY par- 
ticles exceeding 2% are presented in Tables ITTIl 121 The 
complete listing including all decays is available on the 
SPA web-site http : //spa . desy . de/ spa/ 





m, r [GeV] 


decay 13 


decay 13 


h° 


116.0 
4 x 10~ 3 


t~t + 0.077 
bb 0.773 
cc 0.021 


WW 0.067 
gg 0.055 


H° 


425.0 
1.2 


t~t+ 0.076 
bb 0.694 
if 0.052 

f^ff 0.030 


XiX°2 0.038 
X2X2 0.020 
xtXi 0.050 


A 


424.9 
1.6 


t~t + 0.057 
bb 0.521 
tt 0.094 

f£f? 0.036 


X1X2 0.054 
X2X2 0.060 
xtxl 0.163 


H+ 


432.7 
0.9 


v t t + 0.104 
tb 0.672 


xtxi 0.143 
v T Tt 0.071 



Table 9. Higgs masses and branching ratios B > 2% in 
SPSla' from J^f. 



X 


m, r [GeV] 


decay B 


decay B 


X? 


97.7 






X°2 


183.9 
0.083 


e%e* 0.025 
t±t t 0.578 


v e Vc 0.116 
D T v T 0.152 


xi 


400.5 
2.4 


xfW^ 0.582 


XiZ° 0.104 
X2Z 0.224 


xi 


413.9 
2.9 


f^r* 0.033 
i> e u e 0.042 
VtV t 0.042 


XiW* 0.511 
xlZ 0.022 
xlZ° 0.024 
Xih° 0.070 
xih° 0.165 


xt 


183.7 
0.077 


ffv r 0.536 


D t t + 0.185 
y e e + 0.133 


xi 


415.4 
3.1 


e^Ve 0.041 
T^V T 0.046 
fife 0.109 


XiW+ 0.063 
X°W + 0.252 
X+Z° 0.221 
xfh° 0.181 



Table 10. Neutralino and chargino masses, widths and 
branching ratios B > 2% in SPSla' from lb'5t : branching 
ratios for the second generation are the same as for the 
first generation. 
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I 


m, r [GeV] 


decay B 


decay B 




125.3 
0.10 


X?e~ 1.000 




Sl 


189.9 
0.12 


X?e~ 0.925 
X° 2 e~ 0.026 


XiV e 0.049 


Ve 


172.5 
0.12 


X°v e 1.000 




fl 


107.9 
0.016 


Xi r ~ 1.000 




T2 


194.9 
0.18 


Xir" 0.868 
X%T- 0.046 


xr^r 0.086 


V T 


170.5 
0.12 


XiVt 1.000 





Table 11. Slepton masses, widths and branching ratios B > 
2% in SPSla' from \bb^ ; branching ratios for the second 
generation are the same as for the first generation. 



q 


m, r [GeV] 


decay B 


decay B 


UR 


547.2 
1.2 


X°u 0.990 




u L 


564.7 
5.5 


X°u 0.322 


X+d 0.656 


dR 


546.9 
0.3 


X?d 0.990 




d L 


570.1 
5.4 


X 2d 0.316 


X~u 0.625 


ti 


366.5 
1.5 


X°t 0.219 
X°t 0.062 


xtb 0.719 


h 


585.5 
6.3 


X\\t 0.042 
X° 2 t 0.103 


xtb 0.265 
Xtb 0.168 
iiZ° 0.354 
0.059 


bi 


506.3 
4.4 


X?6 0.037 
X_lb 0.295 


0.381 
tiW~ 0.281 


b 2 


545.7 
1.0 


X°b 0.222 
X° 2 b 0.131 
xib 0.028 
xib 0.038 


xr* 0.178 
ttW~ 0.401 


9 


607.1 
5.5 


u R u 0.086 
H L u 0.044 
d R d 0.087 
d L d 0.034 


tit 0.189 
616 0.214 
b 2 b 0.096 



Table 12. Masses, widths and branching ratios B > 2% 
of colored SUSY particles in SPSla' from )65}l: branching 
ratios for the second generation are the same as for the first 
generation. 



(b) LHC and ILC cross sections in SPSla' 

Total cross sections are presented in Figs. |2 - El f° r 
SUSY particle production at the LHC and the ILC. 
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Fig. 2. Total cross sections for squark and gluino pair pro- 
duction at the LHC \27% ££jl for fixed gluino mass (top), 
squark mass (center), and gluino / squark mass ratio (bot- 
tom) [fixed parameters corresponding to SPSla! values]. 
Black circles indicate the SPSla' mass values. The Born 
cross sections (broken lines) are shown for some channels. 
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100 200 300 400 500 600 700 
m [GeV] 

Fig. 3. Generic examples of total cross sections (Drell- 
Yan and Compton production) as a function of the average 
mass for production of sleptons, charginos and neutralinos 
at the LHC \27jf28l . The Born cross sections (broken line) 
are shown for comparison. 
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Fig. 4. Total cross section sections for chargino and neu- 
tralino pair production in e + e~ annihilation \ 73( . The Born 
cross sections (broken lines) are shown for a few channels. 




200 400 600 800 1000 1200 
y/s [GeV] 




200 400 600 800 1000 1200 



yfa [GeV] 

Fig. 5. Total cross sections for smuon and selectron pair 
production in e ± e~ annihilation \74\ - The Bom cross sec- 
tion (broken lines) is shown for comparison. 
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Fig. 6. Total cross sections for fiii pair production in 
e + e~ annihilation for left- and right-handed polarized elec- 
tron (P e - — 0.8) and positron (P e + = Q.Q) beams The 
Bom cross section (broken line) is shown for comparison. 
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